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Real-time controls based on quantum measurements are powerful tools for various quantum proto¬ 
cols. However, their experimental realization have been limited by mode-mismatch between temporal 
mode of quadrature measurement and that heralded by photon detection. Here, we demonstrate 
real-time quadrature measurement of a single-photon wavepacket induced by a photon detection, 
by utilizing continuous temporal-mode-matching between homodyne detection and an exponentially 
rising temporal mode. Single photons in exponentially rising modes are also expected to be useful 
resources for interactions with other quantum systems. 

PACS numbers: 03.67.Lx, 42.50.Dv, 42.50.Ex, 42.65.-k 


Quantum measurement is a basic requirement for 
various quantum protocols. In many of them, real¬ 
time acquisition of the measurement outcomes is ben¬ 
eficial or even required. Real-time measurement en¬ 
ables real-time feedback or feedforward controls of quan¬ 
tum systems. Measurement-based quantum computation 
(MBQC) [U12] is a typical example of such real-time us¬ 
age of the measurement results, where the measurement 
results are fedforward to the next computational step. 
In addition, quantum dynamics conditioned by measure¬ 
ments are often referred to as quantum trajectories or 
quantum filters ma, and form essential parts of devel¬ 
oping quantum control theories. Based on them, basic 
cases of quantum controls, e.g., suppression of noises, 
have been tested with a variety of systems recently [MS]. 

In order to put the quantum measurement-based tech¬ 
nologies further, an important adaptation is to bridge dif¬ 
ferent quantum systems or different quantum variables in 
the measurement system laiiQiiii]. As a general prob¬ 
lem, one may want to connect different systems as a sig¬ 
nal and probe. Aside from the issues of connecting dif¬ 
ferent systems, even when restricting us to the light field, 
there are two fundamental variables reflecting the wave- 
particle duality. One is the continuous variable (CV) field 
quadrature amplitude x = (d + d'*')/\/2, and the other is 
the discrete variable (DV) photon number h = a, where 
d and denote the field annihilation and creation opera¬ 
tors, respectively. CV-DV hybrid architectures have been 
identified to offer crucial advantages, ranging from deter¬ 
ministic and fault-tolerant MBQC [T21IT5] to loophole- 
free Bell inequality test m- Experimental technologies 
for the hybridization are also rapidly developing [17]- 
[20] : however, combining them with real-time quantum- 
optical controls are facing difficulties. 

An obstacle is a mode mismatch between CV-based ho¬ 
modyne detection and DV-based photon detection. To be 


more precise, homodyne detection is sensitive to vacuum 
fluctuation, while photon detection is not. Therefore, 
filtration of the field quadrature, i.e. integration with 
an appropriate weight function /(t) as Y = J f{t)x{t)dt 
is essential in order to remove irrelevant vacuum noises 
and obtain meaningful quadrature information matching 
with the field mode heralded by photon detections. At 
this point, quadrature measurement often fails to be in 
real time. That is critical for quantum feedback or feed¬ 
forward controls since any latency can result in decoher¬ 
ence of the systems. For example, previous experiments 
of CV-DV hybrid MBQC have been limited to Gaussian 
operations [121EO], where the feedforwards are linear, and 
thus filtration can be done after the feedforwards. In or¬ 
der to complete the universal operations of the hybrid 
MBQC [T 2 IIT 5 ] . one non-Gaussian operation is required 
in addition to the Gaussian operations [2TII23] . In non- 
Gaussian operations, the feedforwards become nonlinear 
m, and filtration must be finished before the nonlinear 
feedforwards since they are not commutable. Thus, there 
has been growing demand for the real-time quadrature 
measurement of light fields in DV-based states. 

In this Letter, we demonstrate real-time quadrature 
measurement of a single-photon wavepacket by creating 
the state in an exponentially rising field mode. Quadra¬ 
ture values of the single-photon wavepackets are obtained 
without additional latency in the filtration, thanks to the 
continuous temporal-mode-matching between the homo¬ 
dyne detection and the field mode heralded by a photon 
detection, as explained later. Our experiment is based 
on a specially designed singly-resonant optical parametric 
oscillator (OPO), which we call an asymmetric OPO. Pu¬ 
rity of the created single photon is so high that strongly 
non-Gaussian behavior of the continuous quadrature sig¬ 
nals is obtained in real-time on an oscilloscope without 
any postprocessing of the data or compensation for losses. 
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Single photons in exponentially rising modes are not 
only good for bridging different variables of a light field 
itself, but also for bridging different quantum systems. 
For example, they can perfectly excite two-level atoms 
since they are the time reversals of photons emitted by 
exponentially decaying process of the atoms m- More 
generally, the exponentially rising modes are known to 
be the most efficient modes for interactions with such 
systems as optomechanical oscillators, superconducting 
circuits, etc. [iniiisHii]. Therefore, our single photon 
in an exponentially rising mode can also be utilized as 
single-photon probe fields in quantum trajectories and 
filters [5], or as carriers in quantum networks [28H3Q] . 

There are two possible ways to achieve the weighted 
integration of a field quadrature. One is to tailor the 
envelope of local oscillator (LO) for homodyne detec¬ 
tion to f{t) and then integrate the homodyne signal (op¬ 
tical filtration). The other is to continuously measure 
the field quadrature with continuous-wave (CW) LO and 
electrically perform weighted integration (electrical fil¬ 
tration). The optical filtrations are natural choices for 
pulsed-laser-based experiments [31] |32]. On the other 
hand, the electrical filtrations are used in many successful 
experiments with CW-pumped photon sources j33H36] . 
including some CV-DV hybrid protocols (TSHSO]- The 
electrical filtrations are suitable for such random photon 
sources since they continuously measure the field quadra¬ 
ture. The drawback of electrical filtrations is that they 
have used digital computers for the complicated integra¬ 
tion, which is far from real time. However, when the 
target field mode is an exponentially rising mode around 
time to 

Useit; to) = _ 0 ^ ( 1 ) 

where 0(t) is the Heaviside step function, the filtration 
corresponds to a simple first-order low-pass filter (LPF). 
By adding an LPF to a CW-LO-based homodyne detec¬ 
tion, we can obtain the filtrated quadrature information 

X{tm) = J - t)x{t)dt, ( 2 ) 

which is continuously matched with the exponentially ris¬ 
ing mode. At tm = to, the outcome corresponds to the 
quadrature value of the target field mode /rise(t; to). Such 
a simple LPF can consist of only passive electrical compo¬ 
nents, or a finite-bandwidth homodyne detector itself can 
work as the LPF [37] • Therefore, there is no additional la¬ 
tency in the filtration, and we can overcome the drawback 
of electrical filtrations. The continuous temp oral-mo de¬ 
matching can be extended to more general wavepackets 
as long as they can be matched with simple filters. 

We create a single photon with a heralding scheme, 
where photon pairs are probabilistically generated and 
detection of one photon (idler photon) heralds the other 
photon (signal photon). A conventional doubly-resonant 



FIG. 1. (color online) Single photon generation methods with 
(a) conventional doubly resonant OPO and (b) our asymmet¬ 
ric OPO; Out, outcoupling mirror; APD, avalanche photodi¬ 
ode for idler-photon detection. 


OPO [Fig.j^a)] creates the signal photon in a double ex- 
ponential mode /doubie(i; ^o) = \/^exp(- 7 |t - to\/2), 
where to is the heralding timing and 7 corresponds to 
the bandwidth of the OPO [38]. This is because both 
signal and idler photons are resonant with the OPO, and 
leak from the OPO with the same decay rate 7 . On the 
other hand, our asymmetric OPO [Fig. [^b)] introduces 
asymmetry in the decay rate of the photon pairs by a 
polarization beam splitter (PBS) inside the cavity. Here, 
vertically-polarized (V-pol) signal photon is ejected by 
the PBS, and thus, quickly decays after the detection 
of horizontally-polarized (H-pol) idler photon. As the re¬ 
sult, the wavepacket of the signal photon is approximated 
by /rise(^;^o)- Creating single photons in exponentially 
rising modes have been tested mainly with atomic ensem¬ 
bles [39ti43] . However, these previous experiments have 
reported nearly 50% loss or have been verified just by 
the property of antibunching which is insensitive to the 
loss. Compared with them, the asymmetric OPO has 
several useful features. Transverse modes of the photons 
are cleaned by the OPO so that we can easily get high 
interference visibility with the LO or other beams. Thus, 
we can obtain the information of the state with high pu¬ 
rity. Moreover, our method can be extended to other 
wavelengths by using appropriate nonlinear crystals, and 
7 can be changed by the bandwidth of the OPO. Adjusta¬ 
bility of these parameters suggests that our method can 
create single photons matching with a variety of other 
systems nni na Eg. We also note that we can create 
superposition of Fock states by adding displacement to 
the idler photons [44l[45] . 

The experimental setup is shown in Fig. [^a). We 
use a continuous-wave titanium sapphire laser operat¬ 
ing at wavelength of 860 nm. The OPO, containing a 
periodically-poled KTi 0 P 04 crystal with type-H phase¬ 
matching, is pumped with a 430 nm CW beam which is 
enhances by a build-up cavity. The build-up cavity also 
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FIG. 2. (color online) Experimental setup and the result of 
preliminary experiment, (a) Experimental setup; A, wave¬ 
length; R, reflectivity, (b) Solid trace: estimated temporal 
mode of the single-photon wavepacket. Dashed curve: weight¬ 
ing function of the 3rd-order LPE in the signal line (ii) in 
Fig. I^a). Dotted curve: temporal mode calculated by sub¬ 
stituting the bandwidths of the cavities into Eq. it- 

improves the stability of the transverse mode of created 
signal photon by stabilizing the optical path of the pump 
beam. Since the OPO is pumped with a CW beam, the 
single-photon wavepacket is created at random timing. 
Bandwidth of the OPO is designed to be about 11 MHz of 
half width at half maximum (HWHM). The signal of the 
homodyne detector is recorded with an oscilloscope trig¬ 
gered by the H-pol idler-photon detection at an avalanche 
photodiode (APD). Two filter cavities (FCs) which re¬ 
move photons in irrelevant frequencies are employed in 
the idler line. Their bandwidths are designed to be about 
19 MHz and 36 MHz of HWHM, respectively. Since the 
FCs work as additional Lorentzian filters, the temporal 
mode of the heralded signal photon /Hse(^’^o) becomes 
linear combination of three exponentially rising modes as 

/rtee to) oc El - t), (3) 

n=l 

where 71 , 72,73 correspond to the bandwidths of the cav¬ 
ities, Cl = 1/(72 — 7 i )(73 — 7 i) and 02,03 are its cyclic 
permutations. In this case, the filtration corresponds to 
a 3rd-order LPF. 

For the verification, we compare the conventional dig¬ 
ital postprocessing method and the real-time method by 
simultaneously recording the unfiltered and filtered ho¬ 


modyne signal [signal line (i) and (ii) in Fig. [^a), re¬ 
spectively] . They are recorded during 500 ns around each 
heralding signal with a sampling rate of 2.5 GHz. The un¬ 
filtered signal is analyzed after collecting a series of data, 
using the same manner as [35l|46]. On the other hand, 
the 3rd-order LPF for the real-time measurement con¬ 
sists of only passive electrical components (two chipped 
capacitors and one chipped inductor as in Fig.j^a)), and 
its circuit length is about 5 mm. There is no latency in 
the filtration except for the propagation time in the elec¬ 
tric cables. The event rate of photon detection is about 
1,800 per second, and we collect the data of 18,491 events. 
Bandwidth of the homodyne detector is 100 MHz, which 
is much broader than those of the cavities so that it does 
not distort the detected wavepacket. Since single pho¬ 
tons are phase insensitive, we scanned the LO phase at 
about 100 Hz while recording the data. 

The solid trace in Fig. |^b) is the temporal mode of 
the single-photon wavepacket estimated from the princi¬ 
ple component analysis [35l [46] of the unfiltered homo¬ 
dyne signal. It clearly shows an exponentially rising and 
suddenly falling feature. Fig. [^b) also shows the theo¬ 
retical curve calculated by substituting the bandwidths 
of the cavities into Eq. and the weighting function of 
the 3rd-order LPF in the signal line (ii) in Fig. I^a). The 
mode-match between the single-photon wavepacket and 
these curves are over 99%. Although the performance of 
the real-time quadrature measurement does not depend 
on whether the LPF is Ist-order or 3rd-order, it should 
be noted that we can make the wavepacket closer to sin¬ 
gle exponentially rising mode as Eq. 0 by using more 
broad-band FCs. 

Figure shows the screen captures of the oscillo¬ 
scope displaying the traces of the filtered homodyne sig¬ 
nals without any digital postprocess or compensation for 
losses. Here, the quadrature information correlated with 
the idler field is continuously monitored as the 

real-time signal at tm- The vertical dotted line in Fig. 
indicates the timing of the heralding signals, which corre¬ 
sponds to to in Eq. 0 . Each signal on the white dotted 
line corresponds to the quadrature of the single-photon 
wavepacket. As the number of overlaid traces increases, a 
non-Gaussian dip gradually appears on the white dotted 
line m, which is characteristic of quadrature distribu¬ 
tions of single photons. Continuous change of the distri¬ 
bution with respect to tm in Fig. [^c) reflects the overlap 
between the single-photon wavepacket (t; to) and the 
mode filtrated in real-time Since they are in 

the same shape, the overlap depends on |to — tml- That 
is why we can see the symmetric distribution about to- 
Observing the particle nature of light in real-time by the 
wave-based measurements is in contrast to interference- 
diffraction experiments where wave nature is observed by 
particle-based measurements, e.g., [48]. 

A comparison of quadrature distributions obtained by 
the real-time method and the postprocessing method are 
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FIG. 3. (color online) Screen captures of the oscilloscope dis¬ 
playing the homodyne signal after the LPF mi. They are 
recorded 500 ns around each heralding event, (a) Single trace 
of the signal, (b) Overlaid 32 traces, (c) Overlaid 18,491 
traces, tm is measurement time, and vertical dotted line indi¬ 
cates the heralding signal at time tm = to- Histogram at each 
time is visualized by colors changing from blue to yellow to 
red as the sample density increaces. 

shown in Fig. Figure [^a) is the quadrature distri¬ 
bution obtained by postprocessing of the recorded un¬ 
filtered signal. On the other hand, Fig. Bb) is the 
quadrature distribution of the real-time measurement at 
= ^0 5 which corresponds to the distribution of the 
signals on the dotted line in Fig. [^c). Figure gc) is 
the joint distribution of the quadratures, and it clearly 
shows strong correlation between the quadratures ob¬ 
tained by the two methods. Correlation coefficient cal¬ 
culated from the distribution is over 0.99, which means 
almost the same quadrature values are obtained for each 
event. We also estimate the diagonal elements of their 
density matrices from the quadrature distributions us¬ 
ing the maximum-likelihood method m- The error is 
estimated by the bootstrap method. The single-photon 
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FIG. 4. (color online) Quadrature distributions and estimated 
diagonal elements of density matrices, (a) and (b) Quadrature 
distributions obtained by the postprocessing method and the 
real-time method, respectively. Solid curves are their fitting 
curves, (c) Joint distribution of the quadratures, (d) and (e) 
Diagonal elements of density matrices estimated from (a) and 
(b), respectively. Sectional side view of their Wigner functions 
cutting through the phase-space origin are displayed as insets. 


components of the results of postprocessing method and 
the real-time method are 78.5 ± 0.7 % and 76.8 ± 0.6 %, 
respectively [Figs. El(d) and (e)]. Wigner functions recon¬ 
structed from the density matrices are also shown in the 
insets. They have negative dips at the origin with values 
of —0.181T0.004 and —0.170T0.004, respectively. These 
results clearly show that the highly pure single photon is 
measured with high temporal-mode-matching. 

In summary, we have demonstrated real-time quadra¬ 
ture measurement of a single-photon wavepacket, which 
is created in an exponentially rising mode, with the 
continuous temp oral-mode-matching. We succeeded in 
obtaining the quadrature values of the randomly arriv¬ 
ing single-photon wavepacket with no additional latency. 
From the quadrature distributions, we confirmed the 
high purity of the created single photons, and the high 
temporal-mode-matching between the homodyne detec¬ 
tion and the single-photon wavepacket. Our method is a 
powerful tool for the CV-DV hybrid quantum protocols 
which involves real-time controls. We also expect the 
highly pure single photon in an exponentially rising mode 
will be a useful resource for connecting different systems, 
and will stimulate a variety of studies, from quantum 
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trajectories and filters to quantum networks. 
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